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Abstract 

The problem of the hehx-coil transition of biopolymers in exphcit solvents, like water, with the 
ability for hydrogen bonding with solvent is addressed analytically using a suitably modified ver- 
sion of the Generalized Model of Polypeptide Chains. Besides the regular helix-coil transition, 
an additional coil-helix or reentrant transition is also found at lower temperatures. The reentrant 
transition arises due to competition between polymer-polymer and polymer-water hydrogen bonds. 
The balance between the two types of hydrogen bonding can be shifted to either direction through 
changes not only in temperature, but also by pressure, mechanical force, osmotic stress or other 
external influences. Both polypeptides and polynucleotides are considered within a unified formal- 
ism. Our approach provides an explanation of the experimental difficulty of observing the reentrant 
transition with pressure; and underscores the advantage of pulling experiments for studies of DNA. 
Results are discussed and compared with those reported in a number of recent publications with 
which a significant level of agreement is obtained. 
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I. INTRODUCTION 

The helix-coil transition is a central event in many genetic processes of living matter [l, 2|. 
A number of different approaches to describing the helix-coil transition in biopolymers have 
app eared in the literature; many of them based on spin models j3[jl2|. As shown in Refs. 
10 . the frequently used Zimm-Bragg 13 1 and Lifson-Roig 1^ models are particular cases 
of a more general model based on a Potts-like Hamiltonian. Within these spin models the 
influence of solvent on the helix-coil transition has been the topic of many studies jl^ . This 
is because biopolymers inside cells act in the presence of water with dissolved metal ions 
and other low and high molecular weight compounds Many of the approaches to 

describing solvent interactions are quite diverse and often intractable analytically. 

Solvents differ by the mechanisms of their action on polymer chains dissolved in them. 
For instance, water is a natural solvent for biopolymers and in addition to forming hydro- 
gen bonded networks in the bulk, water also has the ability to form hydrogen bonds with 
biopolymers. Water is the most important biological solvent with very interesting properties 
due to the large number of anomalies present in the water phase diagram that still are a 



matter of debate [16|. 

With a coupled Ising-Potts model Vause, Walker and Goldstein [l3| achieved significant 
success describing lower critical solution points in hydrogen-bonded mixtures. Accounting 
for these studies, during the last few decades slightly different approaches for describing 
the water phase diagram were introduced by Stokely, Debenedetti, Stanley and others jl8 |. 
Within a Bell-Lavis spin model, the possibility of a reentrant phase diagram between low 
and high density phases of water was recently pointed out [l^. 

The key success in these studies was due to the proper accounting for the fact that 
hydrogen bonding between two species takes place only at special orientations of water 
molecules, reflecting the tetrahedral symmetry of water systems. In principle, applications 
of spin models allows capture of the most important features of solvents with the ability for 
H-bond formation, i.e. the directional character of hydrogen bonds and the large entropic 
changes that occur due to this directionality. Especially useful in this sense is the Potts 
model. 

In view of the above-mentioned facts, for description of conformational transitions in 
biopolymers it seems natural to exploit spin models for both the polymer and solvent. 
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Within this context, a detailed description of the water anomahes is often irrelevant for 
two reasons. First, the most interesting and biologically important events in biopolymers 
(helix-coil transition, protein folding) occur under conditions, where water is in the bulk 
liquid state far from freezing or other critical points. Second, it is widely accepted that 
polymer hydration is short-ranged in a direction normal to the polymer axis, and polymer- 
water interactions affect only one (two, in some cases) layer(s) of the water network 20l |21|. 
Conversely, it is very important for consideration of polymers to properly account for the 
possibility of water-polymer hydrogen bond formation. 

Within a spin description Goldstein attempted to include both polymer and simplified 
water models j^. His approach included different Potts variables for descriptions of states 
of repeated units of polymer and solvent molecules. He showed for any ratio of the two 
energies, provided that the energy of polymer-solvent hydrogen bonds is larger than the 
polymer-polymer hydrogen bond energy, the possibility of the helix-coil transition at high 
temperatures, and the reentrant coil-helix transition at low temperatures. 

The helix-coil transition is a constituent part of protein folding and is closely related to 
cold denaturation as well. Currently there is no agreement on a general molecular mech- 
anism that results in both cold and heat denaturation in proteins. Accordingly, there is 
no generally accepted approach or model that allows descriptions of both the direct and 
reentrant transition on the same footing. 

A rather interesting attempt linking the changes in secondary and tertiary structures of 
polymers was reported in Ref. 22|. It was shown that reentrant isotropic- nematic phase 
transitions can be mediated by helix-coil transformations within individual liquid crystal 
molecules. Reentrance of the isotropic phase was shown to be driven by the inverted helix- 
coil transformation. It is well-known that the coil-globule transition in polymers strongly 
depends on the rigidity of the polymer chain j^. Therefore it is not surprising that by 
altering chain rigidity via the helix-coil transition it is possible to tune the coil-globule 
equilibrium. We reach the conclusion that the reentrant helix-coil transition, arising from 
interactions with water through the reentrant rigidity effect, could be the origin of reentrant 
folding and cold denaturation. 

In the framework of a zipper model, Hansen and Bakk (HB) have recently accounted 



for competing effects between protein and water that lead to cold denaturation 



23 



24|. 



Another attempt of describing the influence of water on folding was reported by the authors 
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of Ref. [1^ who included directional features of water hydrogen bonds using Potts variables 
and modeling folding as Go-like. They showed the possibility of cold denaturation depended 
on the ratio of water-water and water-protein hydrogen bond energies. It was concluded that 
effects of pressure on water density are key to understanding cold denaturation in proteins. 
The reported phase diagram agrees quite well with the experimental data obtained on bovine 



pancreatic ribonuclease A. Among the vast arsenal of possibilities, the authors of Ref. 



12 



used one of the simplest models of water 
over solvent degrees of freedom. 



and still were unable to analytically average 



De Los Rios and Caldarelli 



25| used a similar description for water to account for cold 
denaturation in proteins. They borrowed some simplifications from the HB protein folding 
model and concluded that the effective attractive interaction between hydrophobic species 
should depend on temperature. Near the end of their paper they cautioned that care should 
be exercised when defining generally valid effective interaction potentials among amino-acids, 
since during folding the protein exists in an ever changing conformational environment and 
hence strongly depends on the amino-acid distribution and on interactions with water. 

Within the assumption that the transition free energy of a protein, i.e. the free-energy 
difference between the native and denatured states, is a quadratic function of pressure and 



temperature, Hawley et al constructed a melting theory [26|. Based on their analysis the 
resulting bell-shaped 3D curve {AG{P,T)) was sliced at zero surface to obtain an elliptical 
phase diagram of the transition between the native and denatured states. Qui te g ood agree- 



ment was achieved with some of the published experimental results (see 45 



for review 



and comparison with experiments). The same (Hawley) theory was also applied to describe 
the DNA reentrant melting transition j^^] . 

In a rather different context, to produce both cold and heat denaturations, Riccio, As- 
colese and Graziano 



28 



discussed importance of considering the free energy difference be- 



tween the transition states as a quadratic function of temperature. With the same goal, the 



authors of Ref. 



29 



inserted a temperature-dependent hydrophobic attraction as a curve with 
a maximum, without justifying the origin of such dependence. 

Polypeptides and proteins are identical from the viewpoint of their interactions with wa- 
ter. The same interactions must also affect secondary structure formation in DNA. Although 
there are large structural differences between polypeptides and polynucleotides both macro- 
molecules share some important similarities, the most relevant ones being the stabilization 
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of helical structure by hydrogen bonds and having water as the most important (native) 
solvent ^. Generally, when water-DNA interactions are considered the primary focus is on 
screening negative charges of DNA phosphates by water and dissolved counterions. Whereas 
much less attention is usually paid to formation of water-DNA hydrogen bonds. For many 



years the reentrant transition in DNA was elusive and not observed experimentally 30 |. 



More recently experimental evidence for the reentrant melting has been obtainec 



high pressure measurements 



21, 



3lj and single-molecule pulling experiments 



32 



through 
31. Our 



approach provides the theoretical basis from which the reported experimental observations 



of pressure induced denaturation in DNA 



3l| can be interpreted and provides an explana- 



tion of why the reentrant melting transition of DNA is observed in experiments at much 
higher pressure values than is the case of proteins. 

While it is clear that competition between polymer-polymer and water-polymer hydro- 
gen bonds has to be taken into account, the microscopic origin and consequences of such 
competition remains unclear. Analysis of this point stimulated the investigation reported 
here. The aim of this paper is to join current theoretical approaches for characterizing 
conformational transitions of biopolymers in solvent environments where the possibility of 
hydrogen bonding with solvent is explicitly considered. To this aim, results of previous 
studies are built upon and the influence of solvent interactions on the helix-coil transition 



is considered within a particular Potts-like model j9|. We show how orientational degrees of 
freedom of water can be summed out analytically resulting in an effective Hamiltonian term 
with a temperature-dependent interaction energy. This re-normalized temperature depen- 
dent energy is related to the free energy of the melting transition, and the phase diagram for 
the helix-coil transition is reported. Our approach provides a description of temperature, 
pressure, pH, osmotic pressure and denaturant effects on an equal footing, thus providing a 
convenient framework within which to investigate rather complex situations. 



^ Not surprising since both macromolecules are constituent parts of the genome system of a ceU and function 
in the same aqueous environment. For example, the H-bonds between DNA strands are locally broken to 
access the genetic information and the protein is synthesized according to the genetic code read in the 
same environment. 
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II. THE GMPC MODEL WITH SOLVENT 



38 



10 



34 



For this work consider the Generahzed Model of Polypeptide Chains (GMPC) 
with the following Hamiltonian, 

N N 2m. 

-/3i/totai ({7^}, M) = jY.^I^^ + (i - '^^^0 !)• (2.1) 

1=1 1=1 j=l 

The first term, 

N 

- /3i^polymcr {H}) = ^ S^^^ (2.2) 

4 = 1 

is the Potts-like interaction between different parts of the polymer. Here /3 = is inverse 
temperature, the number of repeated units, and J = U /T the temperature-reduced energy 
of hydrogen bonding. The short-hand notation is exploited, eg. 

A-l 
k=0 

where 6{x, 1) stands for the Kronecker symbol and 7; = 1, . . . , Q is the spin variable, which 
can be interpreted as being the number of rotating isomeric states of each repeated unit, with 
values between 1 and Q. The case when 7/ = 1 denotes the helical state whereas all other 
{Q — 1) cases correspond to coil states. Q is the number of conformations of each repeated 
unit and thus describes the conformational variability. The Kronecker delta inside the 
Hamiltonian ensures the energy J emerges only when all A successive neighboring repeated 
units are in the helical conformation. Thus, restrictions on backbone chain conformations 
imposed by hydrogen bond formation are indirectly taken into account. Estimates on the 



structure parameters for polypeptides were A = 3 and Q = 



stranded homo-polymeric DNA, A = 10 — 15 and Q = 3 — 5 



90 (9|,ll0| and for double- 
nj. This GMPC model has 



been shown to be closely related to the Waiko-Sato-Eaton-Munoz (WSEM) model (compare 



the Eq.[10] of ^\ and Eq.[l] of 



or 



40|). The second term, 

N 2m 



i=l j=l 

represents the explicit interaction with the solvent, where / = ^ is the reduced energy of a 
polymer-solvent H-bond. Due to the presence of the term l — S^^^ in Eq. (12. 4p . as opposed to 
the 6^^^^ term in Eq. (l2.2p . the solvent is competing with the polymer for H-bond formation, 



depending on the ratio J/J. As shown below, this competition gives rise to some interesting 
behavior. 

A few remarks are in order regarding descriptions of solvents used here [o], |4l|. They 
are classified according to their relatively high or low molecular weight and according to 
their types of binding, i.e. either reversible or irreversible. Additionally, reversibly binding 
solvents can be divided into major groups according to their mechanism of interaction with 
the biopolymer. That is, those that compete for hydrogen bond formation with repeated 
units of the biopolymer (competing solvent) and those that do not (non-competing solvent). 
For simplicity in the present study only the low molecular weight competing solvent is 
considered. 

It is assumed (i) Polymer- solvent interactions depend on the state (orientation) of solvent 
molecules with respect to the repeated unit, and there are q possible discrete orientations 
of each solvent molecule; (ii) A spin variable /Xj, with values from 1 to g is assigned to each 
repeated unit i. Orientation number 1 is the bonded one, with energy E. 

Some solvents, such as water and urea are able to form hydrogen bonds with nitrogen 



:m nvdrogen 
3ins jll, f H, 



basis of DNA or peptide groups of amino acids in proteins y, la, |42|, |43|. We assume 
repeated units that are not bonded by intra-molecular H-bonds, to be free to form polymer- 
solvent intermolecular bonds. When one intra-molecular H-bond is broken, two binding sites 
for a solvent molecule become vacant. Thus, in the case of polypeptides there are only two 
binding sites per repeated unit, while in the case of DNA there are four (2x2 for A-T pair) 
or six (3x2 for G-C pair) binding sites, so 2m (m = 1,2,3...) spin variables are required 
to describe the interaction between solvent molecules and each repeated unit. The reduced 
energy J of the Hamiltonian in Eq. (12.21) now becomes J = m ^^''''^^'"^ , where Upp and Uss 
are the energies of polymer-polymer and solvent-solvent H-bonds, respectively. 
Using the simple identity, 

exp(j5f)) = l+{e-^-l)St\ (2.5) 
the partition function associated with Eq. (12. ip can be quickly recast in the following form, 
Z = ^ exp (-/3iJtotai({7i}, M)) = 

N q q q 2m (2-6) 

E nti + y^'"'] E E - E n [i + ^(1 - ^1"^)^ ^) 



Here V = e'^ — 1 and R = — 1. 

As detailed in the Appendix, solvent degrees of freedom can then be traced out analyti- 
cally to obtain 

Z{V, R, Q, q, m) = {e' + q- 1)^-^ x Zo{V, Q). (2.7) 

where 



V + 1 = exp[J] = exp[[//T] 



Upp + Us s 

e T 



and where 



a 



1 + 



'^Ups Upp Uss 



1 



1 + 



l + a ^ 



2m 



(2. 



U +U (2.9) 

^ pp \ ^ ss 

is a parameter which sets the balance between polymer-polymer and polymer- solvent at- 
traction. It proved convenient in Eg. (12.81) to introduce the reduced temperature t = 
2T/{Upp + Uss). 

Thus, the partition function of the original model with solvent can be reduced to the 
same model without solvent and renormalized interactions up to a multiplicative prefactor 
that depends on solvent properties only. 



III. RESULTS AND DISCUSSION 

A. Inverse (reentrant) helix-coil transition arising from competition for hydrogen 
bond formation between water and polymer 

Our model incorporates the possibility of an inverse (reentrant) helix-coil transition whose 
origin stems mainly from competition between water and polymer for hydrogen bond forma- 
tion, modeled through different values of a. This is obvious directly from the Hamiltonian 
in Eq. 12. as the two parts it is comprised of contain the same Kronecker 6 symbol with 
opposite signs. As usual the Kronecker symbol is equal to one if an intra-molecular hydrogen 
bond exists, and zero otherwise. 

we note that the o.gi.a, Ha,nnto„,a. ^ »,„ exte,.,o. of the one used Ref g in 
that both an arbitrary range. A, of the intra-polymer interactions and an arbitrary number, 
m, of solvent states have been used, as opposed to A = 3 and m = 2 fixed in Ref y. The 
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same analysis based on the transfer- matrix formalism described in Ref |9| can then be carried 
out. Only final results are reported here. Two quantities are of particular interest. First, 
the helicity degree, the average fraction of hydrogen bonded repeated units, 

1 d 



= {sn 



NdJ 



\nZ- 



(3.1) 



where (. . .) is the usual thermal average over the original Hamiltonian (12. ip . The correlation 
length, corresponding to the spatial scale of correlation damping along the chain at large 
distances, can also be readily computed, 

« - (^) ; (3.2) 

where Ai and A2 are the first and second leading eigenvalues of the transfer matrix (see Refjol 
for details). 
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FIG. 1: (Color online) Dependence of the stability parameter (F + 1)/Q on reduced temperature 
t = 2T/{Upp + Uss) for the polypeptide (PP) parameter set: A = 3, Q = 60 with g = 10, m = 1. 



Curves at different values of a 
solvent-free transition point. 



are shown. Note, the point (1, 1) corresponds to the 



The helix-coil transition occurs whenever C becomes large (see Eg. ( 13.21) ') and hence when 

n n 

Ai ~ A2. On the other hand it was previously shown [8|, |9[ that this corresponds to the 
condition V + 1 ^ Q. Note, in the absence of solvent, with A = 2 this is identical to the 
condition of vanishing AG. This is the change in Gibbs free energy between the helix and 
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the coil state, just as obtained from the Zimm-Bragg model 15| for s = e"^*^/^ = 1, in 
terms of an entropy-enthalpy compensation. 

Figure [1] depicts the stability parameter (V + 1)/Q as a function of the reduced tem- 
perature t = 2T/{Upp + Uss) for several values of the parameter a ranging from a = — 1, 
corresponding to the absence of solvent, to a = 0.3, where the solvent plays a significant 
competitive role. As can be seen from the elementary analysis of Eq. 12.81 as shown in Fig 
1, the {V + 1)/Q curve is monotonic at a values from -1 to 0. At exactly a = the curve 
has a plateau at low temperatures, while at a > a maximum appears. This maximum 
becomes lower with increased a. Therefore in the range of — 1 < a < 0, the curve has only 
one intercept with the f(t) = 1 line at temperatures close to t = 1, indicating a regular 
helix-coil transition. When a > the situation is slightly more complicated, and either 
none or two intercepts exist corresponding to the reentrant transition at low temperature 
and the normal helix-coil transition at high temperatures. This happens because at first, as 
a increases above zero, the water-polymer energy 2Ups becomes slightly larger than Upp + Uss 
and there is a competition between polymer-polymer and polymer-water hydrogen bond for- 
mation. In terms of the Gibbs free energy of the transition, existence of two intercepts of 
the stability parameter {V + 1)/Q with 1, would mean that AG is no longer monotonic with 
temperature, and entropy-enthalpy compensation happens in a more complex, non-linear 
way. In Fig. [T]there exists a well defined value of a ~ 0.04 above which {V + 1)/Q < 1 at all 
temperatures, making the intercept impossible^. This happens because the polymer- solvent 
energy 2Ups becomes so high compared to Upp + Ugs, that mainly intermolecular hydrogen 
bonds are formed and the polymer preferably remains in the coil state. There still can be a 
certain amount of intra-molecular hydrogen bonding. However, as shown below the helicity 
degree is always below one in this regime. 

Eq. 12. 8[ visualised in Fig. [1] represents the result of tracing out the solvent degrees of 
freedom and sets the microscopic basis of the temperature dependence of the transition free 
energy. In the language of the Gibbs free energy, this means that AG is not a quadratic 



function of temperature as typically assumed within the Hawley framework 
Refs. 281 



4J] (also in 



45l-l47l) . Instead, it is a curve with maximum when a > and monotonically 



The limiting value of a is of course different at other Q, q, m values, however, qualitative picture is the 
same. 
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decreases otherwise. 
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FIG. 2: (Color online) Helicity degree and spatial correlation length ^ plotted vs temperature. 
Reduced units were used with the same parameter set as in Fig. [H a = 0.01. 

Additional insight into origins of the transitions can be achieved by considering the he- 
hcity degree, 6*, and the correlation length, ^. As clearly shown in FigjJl at a = 0.01 
is zero at very low temperatures, suddenly increases to 1 remaining at this plateau up to 
temperatures of order 1 (in reduced units), and then drops back to at higher temperatures. 
There are also two peaks of ^ and the temperature values where 9 = 1/2 and peak positions 
of ^ are comparable. These transition points coincide with those determined from the inter- 
cepts with the f{t) = 1 function in Fig. [1] (see the curve corresponding to a = 0.01). These 
quantities unambiguously indicate the presence of a regular helix-to-coil transition at high 
temperatures and an additional reentrant helix-to-coil transition at low temperatures. 

Note that the maxima of correlation lengths for both high and low temperature transitions 

□ 

are equal. This is readily understood in terms of the cooperativity of the GMPC model 
since Q and A are not altered upon a change in a. 

B. Reentrance in polynucleotides 

While the phenomenon of the reentrant protein folding (cold denaturation) transition 
has long been considered in protein melting, much less is known both experimentally and 
theoretically about the reentrant transition in DNA. 
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In a previous study ll| we showed how the model with the Hamiltonian in Eg. I2.2l can be 
apphed to describe the hehx-coil transition in DNA provided that long-range loop formation 
is neglected. Thus, the considerations described above are also applicable to DNA and the 
existence of reentrant denaturation appears rather naturally since hydrogen bonds play an 
important role in DNA as well. This is confirmed by a number of experimental observations 
and seems quite obvious since water is the natural solvent for both. 

Dubins et al reported on pressure effects in double stranded nucleic acid melting [2j 
With the help of Hawley's phenomenological theory 26j they showed that there is a max- 
imum point of the pressure-temperature diagram around 50 °C where the nucleic acid is 
destabilized by pressure at temperatures lower than that and stabilized at higher tem- 
peratures. Using optical absorbance, they reported that the DNA/RNA hybrid duplex, 
poly(dA)poly-(rU) in 20 mM NaCl undergoes a pressure-induced helix-to-coil transition at 
room temperature under elevated pressure. Rayan and Macgregor [3l| reported the spec- 
trophotometric observation of destabilization for poly[d(A-T)] and poly[d(I-C)] at increased 
pressure and various co-solvent concentrations. Thus, there is considerable experimental and 
theoretical evidence for the existence of a reentrant helix-coil transition in polynucleotides. 



C. Phase diagram of the heUx-coil transition and the meaning of parameter a 

As remarked in connection with Fig. [H the helical state is not possible for all values of 
temperature and a. Fig. |3] illustrates this point in the a vs reduced temperature plane, where 
all points above (below) the depicted a{t) curves are in the coil (helix) state. The cases of the 
polypeptide chain (PP), as well as regular DNA heteropolymers of the AT and GC type have 
been considered. Note, the phase diagrams for DNA are much sharper than the polypeptide 
curve. Sharpness of these curves is associated with the higher cooperativity of the helix-coil 
transition of DNA represented by higher values of A (reflecting more rigidity) and by smaller 
Q, corresponding to strong reduction of available conformational space compared to the PP 
case. 

The main obstacle in comparing results presented here with experiments hinges on the 
general difficulty of defining a pressure in spin models. This can however be accounted for 
in the following indirect way. Consider the parameter 

a _ 2Ups - {Upp + Uss) , s 

t~ 2T ' ^'^■'^^ 
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FIG. 3: (Color online) Phase diagrams in variables a and t for the polypeptide (PP) and two 
homopolymeric DNA's: AT and GC. For the PP case the variable set was the same as that in 
Fig. [TJ DNA cases have A = 10, Q = 3, g = 10; m = 2 for AT and m = 3 for GC. The regions 
below the curves correspond to the helix state, while regions above are the coil state. The straight 
line corresponds to the case shown in Fig. [2] and is the a = 0.01 curve. There are two intercepts 
with the PP curve, indicating a low-temperature reentrant transition together with the regular 
helix-to-coil transition at higher temperatures. 




where a = '^^^\t} ^I' It^^'"'^ plays here a role similar to the Flory-Huggins x parameter in 



polymer theory 
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In addition here a is a measure of the relative influence of the 
polymer-polymer attractions {a < 0) versus polymer-solvent attractions (a > 0). 

Usually in order to avoid effects of inter-polymer interactions helix-coil melting experi- 
ments are performed on dilute solutions of biopolymers. Since the volume fraction of polymer 
is extremely small at low polymer concentrations, effects of hydrostatic pressure on polymer 
conformations can only be indirect occuring through changing water properties. It has been 
shown 49| that hydrostatic compression loosens water hydrogen-bond structure near the 



monomer unit. Thus, increasing hydrostatic pressure effectively promotes transfer of wa- 
ter molecules from the hydrogen bonded cluster to non-structured solvent with subsequent 
binding to the macromolecule. Therefore increasing hydrostatic pressure makes binding of 
water molecules with the polymer more thermodynamically favorable than binding to other 
water molecules, resulting in decreased values of Uss- To summarize, it can concluded that 
increased pressure corresponds to increased values of a and vice versa. 

Native, ordered phases are usually found in regions where both pressures and temper- 
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atures are low so that the energies of polymer-polymer and solvent-solvent interactions in 
this region overwhelm the polymer-solvent interaction energies and a < 0. Conversely, at 
low temperature an increase of pressure would result in a disordered coil phase where a > 0. 

This line of logic in view of Fig.[3]suggests the reentrant coil-helix transition we observed is 
the counterpart of high pressure denaturation observed experimentally. In general, stability 
of t he poly mer-water system can be altered in many ways. These include increased pressure 



31 



44- 



a], addition of denaturants 
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5l| , changes in pH [471, l52| , application of osmotic 



pressure 53| or stretching force 32|,|5J]. Since all these factors can alter the balance between 
inter and intra-molecular hydrogen bonds, they can be modeled by changes in a, which in 
its turn mirrors changes of the Flory-Huggins parameters. 



D. Phase diagram of DNA under pulling and unzipping and the inversion of force 
effect at low and high forces 

While there is a long history of pressure, temperature induced melting in biopolymers, 
the role of pulling force as an external parameter is relatively new. Long ago Hawley 
and Macleod |4J] studied effects of pressure on the melting temperature of Clostridium 



perfringens DNA. Unlike the heat induced unfolding temperature of proteins, the melting 
temperature of their DNA did not display any curvature, and was instead a purely linear 
function of pressure over the range studied (0-400 MPa). The first numerical prediction 



of the force- induced reentrant helix-coil transition dates back to 1993 [55|. Since then a 
number of other studies have supported this prediction in DNA. Marenduzzo et al 56[ and 
Orlandini et al js^] predicted existence of unzipping using a computational model. Phase 
diagrams for pulling have been reported by Rudnick and Kuriabova 58| and Rahi et al 59 |. 
Experimental observations of an inversion of force effects on DNA stability (from stabilizing 
at forces < 7 pN, to destabilizing at higher forces) were reported |32| and the corresponding 
phase diagram for pulling was shown in jssl. A similar phase diagram was reported in the 
';heoretical work by Hanke et al |54| in (force, temperature) variables. The phase diagrams in 
33| (experiment, see Fig.5 of jssi) and 54 1 (theory, see Fig. 2a of js^, case where A = 0.01) 



are similar, in that at small forces there are two regions where melting temperature increases 
with force. In these regions force has a stabilizing effect on DNA. 

Our Fig. |3]reproduces the main qualitative features of the phase diagram for DNA melting 
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(see, e.g. Fig. 4 of [33|) and agrees well with Refs. |56l-l59|). However, contrary to those 
studies we did not observe an increase in stability with increasing a on the right side of the 
phase diagram. This might be due to specific features of DNA associated with large loops 



which were omited in our studies 



Another possible source of this discrepancy could be 



the polyelectrolyte nature of polynucleotides. For instance when water-polymer hydrogen 
bonds are broken and water molecules return to the bulk, helix formation is then hindered 
by electrostatic repulsion between negatively charged nearest neighbors that were previously 
screened by water. 



E. Pressure versus force: what is the difference? 



Our analysis provides an explanation as to why it is much more difficult to observe the 
reentrant transition under high pressure compared to experiments using a pulling force. As 
previously mentioned, hydrostatic pressure can only indrectly affect the balance between 
polymer-solvent and polymer-polymer bonds. A very large variation of pressure is required 
to induce conformational changes in dilute solutions of biopolymers. On the other hand, 
applied force operates directly on U^g and t/pp, meaning small changes in pulling force induce 
large changes in a. This because it is much harder to compress several milliliters of water 
enough to produce significant changes in biopolymers occupying a relatively small volume 
of the solution, compared to pulling a single molecule. 

Further going back to the phase diagrams in Fig. [3l in the case of the polypeptide, 
small changes in a shift the reentrant transition point to moderate temperatures making 
experimental observation of cold denaturation possible. In contrast on phase diagrams for 
DNA small increments of a do not result in significant changes of stability and a much 
wider range of experimental pressures is required to induce the transition. As opposed to 
pressure, pulling more strongly affects the hydrogen bonding equilibrium in DNA allowing 
observation of the coil state at experimentally feasible temperatures. 

Together these facts explain why the reentrant transition in DNA is very difficult to 
experimentally observe in pressure-based experiments, but is readily observed in experiments 
using a pulling force. 
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F. Buhot and Halperin approach to helix-coil transition in polypeptides under 
stretching 



Recently, Buhot and Halperin offered an elegant study of extension behavior of helico- 
genic polypeptides. They have explicitly added extension force into the Zimm and Bragg 
description. In one of their papers a step-like behavior of helicity degree versus normalised 



force has been reported [60|, without any model for solvent. A thorough study of Tamas^ 
and Pincus reported no step-like behavoir of helicity degree within the same model 



niro 



61|. 



Later, Buhot and Halperin mention that the appearance of a step on helicity degree is an 



artifact of the approximations they have used 



62l |. The theories by Buhot and 



Ha 



perin and 
, however 



Tamashiro and Pincus are in good agreement with some experimental studies 
there was no intention to study the effect of stretching onto reentrant helix-coil transition, 
and no solvent effects were included into the consideration. 

It could be interesting to experimentally study the pulling effects in the region of tem- 
peratures where the reentrant coil-helix transition takes place. Such a study would justify 
the necessity of solvent model inclusion into the approach offered by Buhot and Halperin. 



IV. CONCLUSIONS 



In this paper within the framework of the Generalized Model of Polypeptide Chains 
(GMPC) the possibility of a reentrant helix-coil transition at low temperatures has been 
investigated. The GMPC is a Potts-like spin model where competition between the tendency 
to form a helix state, as enforced by polymer-polymer interactions J, is balanced by the 
configurational entropy as measured by the remaining Q — 1 coil states. In the presence 
of explicit solvent interactions the equilibrium is modified by adding a competing polymer- 
solvent H-bond formation (strength J), and the solvent configurational entropy measured 
by the remaining q — 1 states. The number of solvent spin variables depends whether a 
polypeptide or a DNA is considered, and in our model is considered an arbitrary variable 
2m in the model, so that both cases can be treated within the same unified model. 

We have shown how solvent degrees of freedom can be traced out exactly to obtain an 
effective GMPC model with renormalized interactions. This is then studied following the 
recipes outlined in past work in the absence of solvent interactions. A low-temperature 
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coil-helix reentrant transiton is found in terms of a parameter akin to the Flory-Huggins 
parameter. This is indicated by the appearence of a sudden drop in the helicity, 6*, and 
confirmed by a second low-temperatures peak in the correlation length, ^. A global phase 
diagram separating helical and coil states has also been illustrated. 

It should be emphasised that the key to our success in qualitatively describing the above- 
mentioned systems and experimental situations stems from the possibility of modeling the 
most relevant solvent features with respect to the helix-coil transition, within a sufficiently 
simple theoretical scheme that allows analytical treatment. 
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Appendix A: Exact integration of the solvent degrees of freedom 



In Eq. (12. 6 P we have 

-^comp.solv. 



2 m 



1 1 



2 m 



R'a - st^) E ^ 1) ^ (/^'' 1) ^ W' 1) + - 

j<k<l 

g2m ^ _ ^(A)^ [2mi?g2"^-i + Cl^R^q^"^-^ + Cl^R^q^'"^-^ + ... + i?^"^] = 



(Al) 



2m 



1 - 5^)) [q + i?]^™ - (1 - 5r)q'- = {q + R) 



i2m 



(A)\ 2m 



\2m 



1 - + 



2m 



where V = e'' — 1 and R = — 1 have been introduced in the main text, and where the 
binomial coefficients = n\/{m\{n — m)!) have been exploited. The last expression can 
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be inserted instead of the solvent part of Eq. f l2.6p . resulting in 



N 



H} i=l 



2m 



2mr(A) 



1 - sr + 



(A) , q Oi 



(q + R) 



2m 



TV 



hi} i=i 



(A)' 



where 



{V + _ exp (m^^^2£+t^)g2m 

Jq + R) 



+ 1 = exp (J) 



2m 



g - 1 + exp (%^) 



2m ' 



Then Eq. fl2.7p is obtained. 
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